Domination of black hole accretion in brane cosmology 
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We consider the evolution of primordial black holes formed during the high energy phase of the 
braneworld scenario. We show that the effect of accretion from the surrounding radiation bath is 
dominant compared to evaporation for such black holes. This feature lasts till the onset of matter 
(or black hole) domination of the total energy density which could occur either in the high energy 
phase or later. We find that the black hole evaporation times could be significantly large even for 
black holes with small initial mass to survive till several cosmologically interesting eras. 

PACS number(s): 98.80.Cq 



The braneworld scenario in which our observable uni- 
verse is a brane embedded in a higher dimensional bulk 
has gained much popularity in recent times. The motiva- 
tion for such a consideration originates from solutions of 
string theory where matter and radiation are confined to 
the 3-brane, whereas gravity propagates also in the bulk. 
The significant consequences of this idea on the physics of 
the early universe is the focus of much current attention 
11. The foremost departure from standard cosmology is 
that there exists a regime during the early stages when 
the expansion rate of the universe is proportional to its 
energy density. The aim is to look for observable impli- 
cations of this behaviour on cosmology. 

Primordial black holes through their evolution and 
evaporation are potentially interesting candidates to- 
wards probing the high energy early stages of brane cos- 
mology. The formation of primordial black holes through 
various mechanisms has been studied in much detail [Q. 
More recently, the diverse cosmological ramifications of 
primordial black holes are being vigorously pursued. Sev- 
eral papers have analysed the implications of a popula- 
tion of primordial black holes in the early universe on 
density perturbations, baryogenesis, reionization and nu- 
cleosynthesis HU. Other investigations have focussed 
on methods of constraining the spectra of primordial 
black holes using observational results on gamma rays, 
background neutrinos, dark matter and quintessence 
Recently, several attempts of describing the formation, 
evaporation and cosmological consesquences of primor- 
dial black holes in the braneworld scenario have also been 
undertaken 

In the present Letter, we consider the evolution of pri- 
mordial black holes formed during the high energy brane 
phase when the expansion rate of the universe is linearly 
proportional to its energy density jl). We restrict our- 
selves to the Randall-Sundrum Type II model J8|. Black 
holes formed during this regime evaporate at a rate pro- 
portional to their effective area times the fourth power 
of their temperature. Using the induced metric on the 



brane, and expressions for the black hole radius and tem- 
perature in terms of the AdS radius I, Guedens et al (?J 
derived an evaporation law which is different from that 
of usual AD black holes. However, for a more complete 
analysis of the evolution of black holes, one must take 
into account the effect of accretion from the surround- 
ing radiation bath as well. This point has been empha- 
sized by Zeldovich and Novikov g] and the ramifications 
of considering accreting black holes have been found to 
be significant in the context of standard cosmology [Q. 
In our analysis, we consider the effect of both accretion 
and evaporation on primordial black holes formed on the 
brane. 

Let us consider the rate of change of mass M of a black 
hole immersed in a radiation bath. The accretion rate is 
proportional to the surface area of the black hole times 
the energy density of radiation. The evaporation rate is 
proportional to the surface area times the fourth power 
of temperature. Thus, M is given by 
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where gbrane represents the effective number of degrees of 
freedom on the brane, tbh and Tbh are the radius and 
temperature of the black hole, tr is the Stefan-Boltzmann 
constant, and pr is the energy density of radiation given 
by 
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where M4 and are the four dimensional Planck mass 
and length, respectively, (t c = 1/2) is the transition time 
to standard cosmology, and the above relations hold for 
t << t c . (Note here that we have neglected the con- 
tribution of evaporation into the bulk given by a term 
proportional to A-KT BH gbuik T%h since it is subdominant 
even for very small black holes M ) . With the above sub- 
stitutions the black hole equation can be written as 
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where (A — 3/(16) 3 7r) and (B ~ 2/ir) are dimensionless 
numberso. The exact solution for Eq.(6) is 
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with M = Mq at the time of BH formation t = to . It can 
be seen that an extremum for the function M(t) exists 
only if (M /M 4 ) 2 < (2A/2B - l)(t /t c ). The relation 
between formation time and mass |7| 



(8) 



shows that such a small initial mass will violate the lower 
bound on horizon mass {M H /M 4 ) > 2 x 10 6 (/ 4 /0 1/3 ob- 
tained from the contribution of gravitational waves to- 
wards CMBR anisotropy ||. This result shows that for 
a radiation dominated high energy phase (t < t c ), ac- 
cretion dominates over evaporation for all times and the 
mass of a primordial black hole continues to grow mono- 
tonically with 




(9) 



Let us now consider the scenario in which a certain 
number density of primordial black holes exchange en- 
ergy with the surrounding radiation by accretion and 
evaporation. Black holes are formed with an initial mass 
spectrum whose range is much debated JTof . In order 
to keep our analysis uncomplicated, the subsequent cal- 
culations are performed by considering an average initial 
mass Mg. We assume that at a time to the fraction of the 
total energy in black holes is /3, and the number density 
of black holes is riBH(to)- Hence, we have 



Pr{to) = Pn{to) + PbhUo) (10) 
Pbh^o) = Ppr{to) = M n BH (t ) (11) 
PR(t ) = (1 - 0)pr(*o) (12) 

The number density of black holes nsjj(t) scales as 
a(t)" 3 , and thus for a radiation dominated evolution 
on the brane, one has (riBH(t)/nBH(to)) = (to/t) 3 ^ 4 
since a(t) oc t 1 / 4 . The net energy in black holes grows 
since accretion dominates over evaporation. The condi- 
tion for the universe to remain radiation dominated (i.e., 
PBuit) < Pn(t)) at any instant t can be derived from 
Eqs.(7,10,ll,12) to be 
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Thus depending on the value of the parameter (3 there 
may or may not be an era of black hole domination in 
the early brane dominated case. We will analyse these 
two cases separately. 

We first consider the situation when the cosmology 
stays radiation dominated up to the time when brane 
effects are important, i.e., t < t c . For this to be the case, 
the initial fraction (3 should be such that 




B+1/4 



(14) 



If the universe remains radiation dominated up to time 
t c , then (for Mq > M4) the mass of a black hole is given 
by Eq.(9). However, the black hole must remain small 
enough, i.e., (M/M 4 ) < (37r/4)(t/t 4 ) for it to obey the 
5.D evaporation law that we have used . These criteria 
can be used to put an upper bound on the initial mass: 
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In the low energy regime pr oc t~ 2 . The black hole 
mass will continue to grow due to accretion up to a cer- 
tain time (say t t ) after which the radiation density be- 
comes too dilute for accretion to be significant. At this 
stage the rate of evaporation is also insignificant since the 
black hole masses have grown by several orders of mag- 
nitude from their initial values. So there ensues an era 
during which the black hole mass stays nearly constant 
over a period of time until evaporation takes over finally. 
This feature is well supported by numerical simulations 
of coupled black hole and Einstein equations for AD black 



1 We assume that the black holes can emit massless particles only and take gbrane = 7.25 as in Q. It is known that the cross 
section for particle absorbtion is increased for relativistic particles ^|, and also for relativistic motion of black holes Q. In the 
present analysis we neglect these effects which tend to increase the value of B by a small amount. 
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holes M . For the present case the accretion rate is smaller 
since the surface area is oc M instead of M 2 for AD black 
holes. Furthermore, the evaporation rate is (ex M _1 ) in- 
stead of M~ 2 . Hence, a critical mass M max is reached 
before evaporation starts dominating from the transition 
time tt onwardsa. Assuming the small mass condition 
((M/M A ) < (37r/4)(*/t 4 )) to be still valid together with 
radiation domination at the time tt, the lifetime for a 
typical black hole can be computed from Eq.(6) (using 
M = M max at t — t t ) and is given by 
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Using Eqs.(8) and (9), the lifetime is given by 
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Guedens et al (7| derived the increase of black hole evapo- 
ration time compared to AD black holes as oc (//ro) 2 , i.e., 
{[t end (M ,5D)}/[t end (M ,AD)} ~ [(t c /t 4 )(M 4 /M )}} 
Our results show that over and above the enhancement 
of black hole lifetime originating from a modified evapo- 
ration law, the black hole lifetime is further increased by 
a factor oc (M max /Mo) purely due to accretion given by 




t end (M, 5D) 
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for t c < t t < t eq . We thus find that the effect of ac- 
cretion in prolonging the life time of black holes is more 
pronounced compared to the 5D evaporation effect for 
small black holes. From Eq.(17) it can be seen that 
for ~ 10 20 , a small number density of primordial 

bD black holes formed with very low or even subplanck- 
ian masees (Mo < M4) could survive up to the era of 
nucleosynthesis and beyond if t t > 10 5 t c . This result 
is likely to have multifarous cosmological consequences 
||[|] which need to be investigated in detail. But for 
black holes with large initial masses (Mo >> M4) which 
form later, there is a lower net mass gain. Hence, for ex- 
ample, if we take t t ~ 10 5 i c , and (l/U) ~ 10 20 , black 
holes formed with M = 10 15 M 4 ~ 10 10 g, will have 
M max ~ 10 15 g, and will be evaporating now. However, 
if the AdS radius is the maximum allowed by the current 
experimental bounds |0|, i.e., (Z//4) — 10 30 black holes 
with initial masses as low as Mq = 10 8 M4 ~ 10 3 g survive 
up to the present era even for tt ~ t c . 



We next consider the case when the initial energy 
fraction (3 is such that the the accretion process causes 
PBH(t) to exceed pa(t) at a time t < t c . If (3 exceeds 
the bound (we still assume « 0(1)) given in Eq.(13), 
then the cosmology enters a matter (black hole) domi- 
nated regime in the high energy phase when H 2 oc p 2 . 
The onset of such a matter dominated era is -Labelled by 
theq which can be written from Eqs.(9-12) ascl 
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of a 



black hole at t} leq is given by 



The mass 
M(t heq /M ) = r 

For t > theq the Hubble expansion is essentially driven 
by the black holes (p = and H oc pbh)- The number 
density of black holes scales as matter (ngjj(t) oc a~ 3 ), 
and thus the scale factor grows as a ~ t 1 ^ 3 . During this 
era, the radiation density pa is governed by the equation 
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where the r.h.s (contribution from accreting black holes) 
cannot be neglected in comparison with the redshiting 
term (pn ~ a~ 4 ). Taking into account the radiation 
dominated expansion for t < theq, and black hole domi- 
nated expansion for t > theq, one has 



n B H(t) 
riBH(to) 



Using Eqs.(9-12) and (21) in (20) gives 




(21) 
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Using the condition of matter-radiation equality at theq, 

i-e., [pB,{theq) = PBH (theq) = M(theq)nBH(theq)], OllC 

can solve Eq.(22) to obtain 
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where we have neglected a term of higher order in theq/t. 
The black hole equation for t > theq is given by 



2 The exact determination of tt requires the numerical integration of the coupled Einstein and black hole equations (see, for 
example M) and is beyond the scope of the present paper. 

' s Although the formation time to of individual black holes is spread out in accordance with the initial mass spectra, we have 
used to to signify the mean formation time in the similar sense of using Mo as the mean initial mass. 



3 



M = Bj~ 
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where C = (3B/ (B + 1/3) Note that we have ignored the 
evaporation term at this stage since it is even more negli- 
gible compared to accretion because of the resultant mass 
gain up to theq- Using M(th eq /M ) — -f B , the solution 
of Eq.(24) can be obtained. Beyond theq the black holes 
continue to gain mass until the two terms in Eq.(24) be- 
come comparable. At this stage one enters the transition 
regime (tt) after which evaporation takes over. Setting 
the r.h.s of Eq.(24) to at t = tt, one obtains 



tt 
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Thus the accretion regime lasts for a brief duration in the 
matter dominated phase, and the maximum black hole 
mass M(t t ) is given from the solution of Eq.(24) by 
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In the evaporating regime (t > t t ), the rate of change 
of black hole mass is given by the first term in Eq.(6). 
This can be integrated to give 
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The energy density in radiation pii(t) whose time evolu- 
tion is governed by Eq.(20) now has a positive contribu- 
tion from black hole evaporation. Making the appropri- 
ate substitutions, one can again integrate pii(t) to obtain 



PR {t) « 3/3p(i ) 7 - 3/4 
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The universe gets reheated as pn(t) increases with time. 
The stage of black hole domination lasts up to a time t r 
when pn(t r ) — nBH(t r )M(t r ). Subsequently, radiation 
domination takes over again. From Eqs.(27) and (28), 



t c _ 3 A M 4 
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domination take place before t c . Requiring t r < t c , one 
gets a lower bound on (3 from Eq.(29), i.e., 
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The evaporation time of black holes in this scenario 
can be calculated from Eq.(27). The mass gain between 
theq and U given by Eq.(26) is negligible if Eq.(30) is 
satisfied and also << 0(1). The black hole lifetime is 
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The contribution of accretion in prolonging the lifetime of 
black holes compared to the modified evaporation effect 
is given in this scenario by 



t end (M, 55) 
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This enhancement is independent of the value of the ini- 
tial mass. Thus in this case the lifetime of the whole mass 
spectrum gets prolonged due to accretion by the same 
scaling factor 7 2S , as distinct from the scenario of radi- 
ation domination throughout where the low mass black 
holes are affected most. Thus, for instance, 1H/I4 ~ 10 20 
and = 1CP 3 , black holes with Mn ~ 10 12 g evaporate 
during the present era. 

To summarize, we have studied the evolution of pri- 
mordial black holes formed during the radiation domi- 
nated high energy phase on the brane. These black holes 
obey a modified evaporation law Q]. We have shown 
that the accretion of surrounding radiation completely 
dominates evaporation as long as radiation domination 
persists. This results in the net growth of mass of the 
black holes. Compared to the case of black holes formed 
with the same initial mass in standard cosmology, the 
black holes in the braneworld scenario evaporate much 
later. We have found that the effect of accretion in en- 
hancement of black hole lifetime could be very signifi- 
cant for even low mass black holes for a wide range of 
parameters. It needs to be emphasized that the effect of 
accretion and its consequences of black holes surviving 
through key epochs in cosmology |||J need to be studied 
in more details. Such investigations could constrain the 
initial energy fraction in black holes vis-a-vis the size of 
the extra dimension. 



The standard low energy cosmology should be radia- 
tion dominated with sufficient radiation temperature for 
nucleosynthesis to occur. Thus the era of black hole 
domination should be over clearly before nucleosynthe- 
sis. Since it is possible for t c jt 4 to be as large as 10 30 
pi, let us simply demand that the onset of radiation 
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